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§1. Introduction 

There are many systems that electron-electron interactions play important roles in the physical 
properties of them. A system of metal with dilute magnetic impurities is a good example. Various 
anomalous phenomena are observed at low temperatures, related to the Kondo effect, which is 
caused by competition between quantum coherence effects and electron-electron interactions.S) 

In the phenomenon of the electron tunneling through a quantum dot between leads, the effects of 
the electron-electron interactions in the dot have been investigated extensively.a'U'u'B 
The importance of the Kondo effect in a quantum dot has been pointed out by some theoreti- 
cians .a'U'Q'tj' The 'Kondo effect in a quantum dot' will give the following features. In low 
temperatures than the Kondo temperature of the system, Tk, the conductance will show the co- 
herent resonant tunneling caused by the Kondo resonance peak in the quantum dot. When the 
temperature rises up than Tk, the Kondo resonance peak vanishes and the local spin freedom ap- 
pears, and so the electron tunneling with inelastic spin scattering process will dominate. Therefore 
the coherent resonant tunneling process will be broken. 

The interference phenomenon is the good indication for the coherency of the electronic state. The 
Aharonov-Bohm (AB) circuits with quantum dots have been investigated to study the coherency 
on quantum dotsJll'0Ei'lll0Ei) Akera had studied the system in which two quantum dots are 
included in AB circuit, and shown that the amplitude of the AB oscillations depends on the 
spin state of the dots.© However, his investigation was restricted only to the high temperature 
region than Tk. Main purpose of the present work is to study how the Kondo effect affects the AB 
oscillations. We calculate the temperature dependence of the AB oscillations in a wide temperature 
range including the crossover temperature Tk. 

In this paper, we consider the system shown in Fig. 1, which is the same system that Akera 
studied. The essential physics of this system can be given by the two impurity Anderson Hamilto- 
nian. Because our interests are related to the change of the coherency of the electronic state, we 
have to calculate the conductance without assumptions on tunneling processes. In other words, 
the assumption, such as the restriction to the 'coherent resonant tunneling process' or to the 'se- 
quential tunneling process', should not be used. Numerical Renormalization Group (NRG) method 
has been known as a useful method to study the many body systems written by the Anderson 
Hamiltonian This method has reliability in the temperature range near and below Tk. In 

this paper, we develop this method to calculate the conductance without assumptions on tunneling 
processes. 

Main results of this paper are as follows: When the temperature rises up, the amplitude of the 
AB oscillations becomes smaller reflecting breaking of the coherency by spin inelastic scattering; 
Various types of AB oscillations appear due to many body effects when the depth of dot's potential 
is changed; AB oscillations have strong higher harmonics components as a function of the magnetic 
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Fig. 1. The geometry of the two quantum dots and leads. Electron tunneling from the one lead to the other lead 
can occur only via quantum dots. The magnetic flux $ is applied between leads and dots to cause AB interference 
effect. 



flux reflecting the crossover of the spin state. 

In section ^, the model Hamiltonian of the system is given and is transformed for numerical 
calculation. The conductance formula is derived based on the linear response theory. At the same 
time, for comparisons, the conductance formula based on the coherent resonant tunneling process 
is derived by using the Friedel sum rule. In section ^, numerical results of the conductance are 
shown. In section summary and discussion of this study are given. 

§2. Formulation 

2.1 Model and transformation 

For simplicity, we consider following situation for our system in Fig. 1. Energy separations 
between states of the single-electron orbits in a dot are much larger than the other energy scales: 
the Coulomb repulsion energy [/, temperature T, and the level width A. This situation corresponds 
to the assumption that only one single-electron orbit in the dot contributes to the current. Two 
electrons can occupy in this orbit because of the spin degeneracy, but the second electron need 
excess energy U due to the Coulomb interaction between two electrons. The orbit in the dot 
couples to those in leads, and this causes the electron tunneling. Tunneling from the one lead to 
the other lead can occur only via quantum dots. The essential properties of the system can be 
written by the Anderson model as follows; 

H = H, + Ha + i7i_d, (2.1) 

^1 = E ^afc<x«L°fc- + E (2.2) 
iter PO" 

Hd= ^d,idl^dia + Y Uidl^di^dl^dii, (2.3) 

i=l,2 a i=l,2 
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ka 

pa 

(2.4) 

where a^^ is the annihilation operator of the electron with spin a and the state k in the left lead A, 
bpcr is that of in the right lead B, and di^a is that of on the orbit in the i-th dot, respectively. The 
quantity £d,« is the single-electron level of orbit in the i-th. dot, and is the effective depth of dot's 
potential. The quantity Ui is the Coulomb repulsion energy in the i-th dot. The terms Hi and 
denote the leads and the dots, respectively, and -ffi_d is the electron tunneling term between leads 
and dots. The factor e 4 is due to the magnetic flux, (ip = 27r$/$0) where $ is the magnetic flux 
enclosed by the circuit, and $0 is the magnetic flux quantum, hc/e.) 

We transform the model Hamiltonian to fit numerical calculation. We consider only the situation 
that £(1,1 = £d,2 = £d Sknd Ui = U2 = U for simplicity. Furthermore, we assume that the geometry 
of the leads has symmetry with respect to the interchange of dot 1 and dot 2. Moreover we also 
assume the symmetry with respect to the interchange of lead A and lead B. Electron motion in 
the leads are restricted only to the onc-dimcnsional motion along the lead. (The state index k 
is reduced to k.) Neglecting the k dependence of the tunneling matrix, we put the relations, 
'^ak,i = "f^bp,! = i'aA;,2 = 'y6p,2 = V. Finally, each term is rewritten as 



Hi = J2 £fc(as,fca«s,fc<x + al k^a^M)' (2-5) 

ka 



+ jK,T + "ea)K,T + "oa) 

{^)aia2 ' ('^)CT3(T4^e,<Ti'^e,CT2^o, 0-3 '^0,0-4 

<Tl(T2Cr3<T4 

--^i4,r4,ido,ido,r + dl^dl,i'^c,idc,r), (2.6) 
H\-d = 2 ^(y cos ^dl^as^ka + V sin ^dl^^as,,ka 

ka 

+h.c.), (2.7) 

where a^^ka = [o-ka + bka)/V^ is the symmetric combination of the lead orbits, and aa.fco- = {o,ka ~ 
bko-)/V^ is the anti-symmetric combination of those, and d^^a = {dia + d2a)/V^ , c?o,o- = —^{dia — 
d2a)/V^ are those of the dot orbits, respectively. The quantity a is the Pauli matrix. The problem 
is reduced to solve the two-channel Anderson model as seen in eqs. (2.5), (2.6) and (2.7). Hereafter, 
we call the orbit denoted by dea as 'even-orbit' and the channel of it as 'even-channel', and those 
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of do o- as 'odd-orbit' and 'odd-channel'. We denote them by / as / = e or / = o. 

2.2 Conductance formula 

In this paper, we restrict ourselves to the linear response conductance for the applied bias voltage. 
We derive the conductance formula without assumptions on electron tunneling processes. We 
also derive the conductance formula with assumption of the coherent resonant tunneling process. 
Calculated results are compared at very low temperature. 
2.2.1 Conductance formula at finite temperature 

Now, we must calculate the conductance without using approximations for tunneling processes 
such as the sequential tunneling process picture or the coherent tunneling process picture. We 
define an electric current from lead A to lead B as follows; 

I^-edMpI^, (2.8) 

where — e is the charge of an electron, and (A^a) is the expectation value of the time differentiation 
of Na = J2ka Ofco-Ofco-, the electron number operator in the lead A. The quantity (Nb) is that of the 
lead B. We consider the situation applying the voltage 2V between lead A and lead B. This situation 
is given by adding the term H' = N^eV — N-qcV to our Hamiltonian H, eq. (2.1). Assuming that 
the voltage difference 2V is small, and using linear response theory derived by Kubo,@) we obtain 
the following expression for the conductance formula (see Appendix); 



with 



4 Z ' ' 



X 



n,m 

2 



{n |Aa - A^B m) 
x5 (u; - {En - Em)) , (2.10) 

where Z = J2n^~^^" the partition function, (3 = 1/T, T is the temperature of the system. 

The quantities (Aa) and (Ab) can be expressed by localized operators near the dots as shown 
in Appendix. So, this expression is suitable for the numerical calculation by the NRG method, 
though it needs the delicate limiting process u; — > 0. 

2.2.2 Conductance formula at zero temperature based on the coherent resonant tunneling process 
Here we derive the conductance formula at zero temperature for comparison with numerical 
results from eqs. (2.9) and (2.10). The conductance at zero temperature is given by 

Gf = — ^ l^fin,int(eF)|^ 
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xp\e^)Ws,^,, (2.11) 

where initial and final states of transition are denoted as int = {ka, Sjnt) and fin = {pa , 5fin). The 
quantity k denotes the electron in the lead A, and p denotes that in the lead B. By S\rit and S^^ we 
represent initial and final states on the quantum dots, and Ws-^^^^ is the probability of state in Si^t- 
The quantity T^n,\nt is the transition matrix, and p{e-p) is the density of states on the Fermi energy 
in the leads. The transition matrix is written by using Green's function Q asT = V + VQV, where 
V = -ffi-d in this case. 

At zero temperature, the system written in eq. (2.1) is expected to be in the local Fermi liquid 
state. Electron tunneling processes through the quantum dots will be given by only the coherent 
resonant tunneling process, i.e., ^int = S'fin- The conductance Gp is written as follows; 

2g2 

Gf = — lA^GeA^F + iO) - AoGoA^F + iO)|^ 
h 

(2.12) 

with 

Ae = 47r|l^|V(£F)cos2|, (2.13) 
Ao = 47r|y|V(£F)sin2|. (2.14) 

In the local Fermi liquid state, Green's function Gi^ai^F + iO) in the l-th channel (/ = e, o) satisfies 
the Friedel sum rule as followJ^l'li); 

GiA^F + iO) = 7 ^ 7, (2.15) 

where Si is the phase shift for the l-th channel. Using the relation between the phase shift 5i and 
the ground-state occupation number {tlq^i) on the l-th orbit, 

(^0,/) = 2-, (2.16) 

TT 



2e2 . 2 r^. 



we get the conductance as follows; 

Gf = ^ sin^ ||((™o,e) - (no,o»| • (2.17) 

Although the expression (2.17) seems to be different from eq. (2.9), it can be derived from eq. (2.9) 
when the system is in the local Fermi liquid state and T = 0. 

The expression (2.17) is contrast with the single dot case. In such a case, the conductance is given 
by (2e^//i) sin^(7r(no)/2), and it has almost the maximum value 2e^/h when there is the Kondo 
resonance peak on the Fermi energy, i.e., (no) ~ l.i'i'& On the other hand, the conductance has 
very small value in the present case as seen from (2.17) if the even and odd components have 
simultaneously the Kondo resonance on the Fermi energy, i.e., (?io,e) ~ {^o,o) ~ 1- This is caused 
by the interference cancellation between processes through the even and odd orbit states. 
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2.3 Excitation spectra 

We will calculate excitation spectra other than P" (co) to get insights of the electronic states of 
the dots system. The single particle excitation spectrum for the l-th orbit {I = e, o) is defined as 
follows; 



ftH = ill 



z 

X (\{n\dl\m) ^ 6{uj- {En- Em)) 



+ \{n\di\m)\' 6 {io + {En - E^))) . (2.18) 
The magnetic excitation spectrum at T = is calculated as follows; 

X:,e(^)=EEK^I(^e,. + 5o,.)|Gr)|2 
n Gr 

x6{u;-{En-EGr)), (2.19) 

where Si^z = {d]i^di^^ — d] ^di^i)/2 is the spin operator on the ^-th orbit and Gr denotes the ground 
state of the system. The energy of the peak position will reflect the characteristic energy of the 
spin fluctuation. In this paper, the energy of the peak position of x'mei^) used as the definition 
of the Kondo temperature, Tk. For the four fold degenerate model, this definition gives rather well 
approximate value of Tk defined by the usual definition based on the susceptibility® 

§3. Numerical Results 

We have several parameters for the model. The band width of the leads, D, is chosen to be 
energy units, D = 1. The parameters e^, U,A = 4:Tt\V\'^ p{ep) and T are given in units of D. 

In this paper, we fix the Coulomb repulsion energy to be C/ = 0.10 and the tunneling intensity 
A = O.OSvr. It is enough to study the only > —0.5U cases because the conductance as a 
function of Ed is symmetric with respect to Ed = —0.5U. Then the energy level Sd is varied 
from —0.05 = — 0.50C/, corresponding to the electron-hole symmetric case, to 0.1 = l.OU. The 
conductance is the periodic function of the magnetic fiux cp and moreover it is even function of 
if for the present model, so we calculate the conductance only in a half period, < < vr. 
Temperature T is varied from 5.8 x 10^^ to 3.1 x 10^^. (See Table I and II.) 

Before going to detailed discussions of the numerical results, we stress the following fact. The 
Kondo effects in different channels seem to occur as if they are independent though the original 
Hamiltonian includes interaction terms between electrons in even and odd orbits. (See eq. (2.6).) 
As shown later in § |3.2| , the Kondo resonance peaks of the single particle excitation in even-channel 
do not have structure reflecting the energy scale of the odd-channel Kondo effect, and vice versa. 
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Table I. 


Parameters ea, U 


, A s 47T\Vfp{eF) in 


units of the band width D. 


£d 




U 


A/tv 


-0.050 = -0.50U 
-0.030 = -0.30U 
-0.025 = -0.25U 
-0.020 = -0.20?7 

+0.000 = +omu 

+0.050 = +0.50U 
+0.100 = +1.00[/ 




0.10 


0.03 


Table II. Temperatures denoted as T, 


in this paper. These 


are given in units of the band width D. 


Ti T2 




T4 


Ts Te Tr 


5.8 X 10"'^ 1.7 X 10"^ 


5.2 X 10"^ 


1.6 X 10"^ 


4.7 X 10"® 1.4 X 10"* 4.2 x 10"* 




Ts Tg 


Tio 


Til 


Tl2 Tl3 


1.3 X 10"^ 3.8 X 10"* 


1.1 X 10"^ 


3.4 X 10"^ 


1.1 X 10"^ 3.1 X 10"^ 



3.1 Conductance at zero temperature Gf{v) 

The conductance at zero temperature is written by using the difference of the occupation numbers, 
{i^o,e) — {no,o), as shown in eq. (2.17). The expectation value of the occupation number (rzo,/) is 
calculated from the wave function obtained by the NRG method, and is shown in Fig. 2. In the 
electron-hole symmetric case, £d = —0.050, the expectation values satisfy the relation (no,e) = 
{i^o,o) = 1-0 for any (p. When moves away from = —0.050, both (no,e) and (no,o) deviate 
from 1.0 and they depend on the flux tp. We note that the hybridization strength for the odd- 
orbit, Aq = 47r|yp/9(£d) sin^((/p/4), is very small in the ~ region. Therefore the small change 
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Fig. 2. The expectation values of the occupation number of the even-orbit, (no.e), and the odd-orbit, (no.o), in the 
ground state as a function of <^ for various cases of ea- The flux is given in units of $o = ftc/e. The value (no,o) 
at 1^ ~ decreases rapidly as sa increases between ea = —0.025 to —0.020. 



of the parameters in a critical region causes drastic change of the ground state properties of the 
odd-channel. Reflecting these, (no,o) at ~ decreases rapidly as ej increases from —0.025 to 
—0.020. At the same time the Kondo resonance peak in the odd-channel rapidly moves away from 
the Fermi energy. 

We have two regimes for the (^-dependence of (no,o)- In the ea < —0.025 case, (no,o) has value 
almost 1.0 at (/? ~ 0. It initially decreases with increasing and then increases gradually. In the 

> —0.020 case, (ng^o) has very small value at ip 0. It increases with increasing (p. In contrast 
to the behavior of (no,o), the occupation number (no,e) shows mild dependence on and (p because 
the hybridization strength Ac is relatively large in < < vr region. 

In Fig. 3, we present the numerical results of the conductance at zero temperature Gf(<^) 
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for various ea- The AB oscillations vanish completely in the electron- hole symmetric case,c3) 



Fig. 3. The conductance at zero temperature, Gp(<^), for various cases of Ed- The units of Gf(v) is 2e^ jh. In 
the electron- hole symmetric case (ea = — 0.50C/ = —0.050), AB oscillations vanish completely. As £d shifts from 
£d = —0.050, the AB oscillations appear. The conductance Gp(^) near c/p = increases rapidly when £d increases 
from £d = -0.025 to -0.020. 



ffd = —0.050, because (no,e) = (^^0,0) = 1 as noted in Fig. 2. When moves away from = —0.050, 
the AB oscillations appear because the quantity (no,e) — ("-0,0) has non zero value and depends on 
the flux Lp. 

We note that Gf(v5) near ~ increases rapidly as increases between = —0.025 to —0.020. 
This phenomenon is reflecting the behavior of (no,o)- In the region of (n.0,0) ~ 0, the Kondo 
resonance peak of the odd-channel moves away from the Fermi energy. Therefore, the contribution 
of the resonant tunneling via the Kondo resonance of the odd-channel becomes very small. From 
this we can see that the conductance Gp{ip = 0) take value near 2e^//i when only the coherent 
resonant tunneling of the even-channel contributes. 

To conclude, the AB oscillations of the conductance at zero temperature show two types of 93 
dependence reflecting the change of the ground state properties of the system. In the Ed < —0.025 
case, Gf(v') is very small at (/? ~ 0. It decreases initially and soon increases with ip, and then it 
decreases near c^/vr ~ 1 after showing maximum. In other word, the AB oscillations have higher 
harmonics components. On the other hand in > —0.020 case, Gf{ip) has large value at <^ = 0, 
and it decreases with ip. 
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3.2 Conductance at finite temperature G{(p) 

In this subsection we show the numerical results of the conductance G{tp) at various temperatures 
calculated by eqs. (2.9) and (2.10). 

In Fig. 4, we show examples of the spectra, P" {uj)/uj, calculated by the NRG method. The 



Fig. 4. Examples of the spectrum, P calculated by using the NRG method. The parameters are = —0.025 

and T = T^i^ 4.7 x 10"^). The symbols give the NRG raw data points, and lines are given by spline interpolation. 
The edge of the lines give the lowest energy data points for the interpolation, |a;niin| = 4.0r5. The data points 
u/T < a ~ 1 should be discarded in the NRG calculation, see the text. 



calculated raw data are plotted by symbols, and the lines give the values averaged by the spline 
interpolation method. The data points with energies uj /T < a ^ 1 are expected to be erroneous 
because they are calculated from transitions with excitation energy (T ~ \tL\) larger than w, where 
I til is the hopping matrix characterizing the low energy scale of the NRG calculation.0ill' In the 
interpolation, these low energy data are discarded. The value at Wmin = 4.0r is substituted for 
the limiting value of w — > 0. As shown from Fig. 4, this process does not lead the ambiguity of 
the limiting value of P"{uj)/uj for the low temperature cases T <^ Tk. But it causes ambiguity 
for the T > a' ~ Tk cases. The quantitative accuracy of the numerical value of the present work 
should not be so trusted for higher temperature cases. However, it seems not bad for the qualitative 
discussions. 

In the following subsections, we show two typical cases for the numerical results of the conduc- 
tance at finite temperature. 
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3.2.1 Deep case (e^ = -0.025; 

Numerical results of the conductance G{ip) at = 
in Fig. 5. 



—0.025 for various temperatures are shown 



Fig. 5. The conductance G{ip) at Sd = —0.025 for various temperatures. The regions given by the broken hues are 
the magnetic state and the regions given by the sohd hues are the non-magnetic state. For the definition, see the 
text. The curves of the conductance given by the solid Une region for Ti < T < Tj cases coincide to Gf(¥') of 
Ed = —0.025 case in Fig. 3. TunneUng processes are expected to be dominated by the coherent resonant tunneUng 
process on soUd Unes. On broken Unes, contributions from the process with the spin inelastic scattering will be not 
so small. The conductance rapidly changes at i^/tt ~ 0.10 in the temperature cases Ti < T < Tj. For the origin of 
this phenomena, see also the text. 



First, we discuss the results in the temperature range Ti < T < T-j. The conductance G{ip) 
changes rapidly as if changes between 0.1 < tp/iT < 0.3. In addition, it has value near 2e^/h in 
ip/ir < 0.1 region contrasted to the result of Gf{ip) for £(j = —0.025 in Fig. 3. In the region of 
(p/ir > 0.3, the curves of the conductance G{(p) overlap to each other, and also coincide with the 
result from G-p{ip). 

To get insights the origin of this behavior of the conductance, we investigate the electronic state 
of the system. The one-particle excitation spectra 7rAePe(a;) and 7rAoPo(u^), and also the current 
spectra P" {ui) for various flux cases at T = Ti are shown in Fig. 6. The same quantities at T = 
are shown in Fig. 7. At T = Ti, we have the Kondo resonance peaks in the even-channel for all flux 
(p cases. (Notice that the abscissa have a logarithmic scale.) On the other hand in (^/tt < 0.1 cases 
for the odd-channel, the widths of the Kondo resonance peaks are small compared to the lowest 
energy calculated in the Fig. 6. The width gradually increases as flux (p increases from ip/ir ~ 0.10, 
and the complete Kondo resonance peak appears in the ip/ir > 0.15 cases. 
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Fig. 6. Excitation spectra 7rAePe(w), 7rAoPo(u;), and P ai T = Ti and Ed = —0.025. The Kondo resonance 

peaks on Fermi energy always exist in Peico)- On the other hand, There are no Kondo resonance peaks in po(w) 
when (p/tt < 0.1. The Kondo resonance peaks in po(w) gradually grow up as flux increases from i^/tt ~ 0.1. The 
peaks appear completely in the ip/n > 0.15 region. 



The evolution of the Kondo resonance peak in a certain channel means that the electronic state 
in the channel is in the spin singlet state. In this case, the excitation properties will be described 
as the Fermi liquid state. Therefore, the electron tunneling processes through that channel will be 
given by only the coherent resonant tunneling process. 

The even-channel state is always in the Fermi liquid state at T = Ti for whole range of ip. The 
electronic state of the odd-channel falls into the Fermi liquid state for (p/ir > 0.15 at T = Ti. On 
the other hand for the cases (/j/tt < 0.1, the spin state of the odd-channel at T = Ti is expected 
to be magnetic. Fig. 7 shows the spectra at T = Ts. The intensity of the Kondo resonance of 
odd-channel with ip/ir = 0.141 decreases remarkably from the lines in Fig. 6. At the same time 
the conductance increases. 

It may be advisable to consider the Kondo temperature, Tk, for given case. Now, we define 
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Fig. 7. Excitation spectra 7rAcpe(<i'), 7rAoPo(w), and P {u))/ijiJ at T = T3 and Ed = —0.025. The Kondo resonance 
peaks in po(u') appear in the Lp/ir > 0.25 region. The region of ifi that the odd-channel being magnetic state becomes 
wider when the temperature rises as seen from comparison with Fig. 6. 



Tk(v?) as an energy of the peak position of the magnetic excitation spectrum at zero temperature 
Xm,e(^) given by eq. (2.19). The result is shown in Fig. 8. In the figure we also show the 
spectrum Xm,e(^) several values of ip. The energy of the peak position is expected to reflect the 
characteristic energy of the spin fluctuation on the odd-channel. For ipj-K = 0.125, we have peak 
at a;/£) = 5.2 X 10~^ and the energy of the peak position increases rapidly as (p increases. For the 
cases ipj-K < 0.125, the peak position shift to very low energy, and thus we can not show it in the 
frame of the figure. 

The judgments of the boarder between the magnetic state (broken lines) to non-magnetic state 
(solid lines) in Fig. 5 were given by using Tk((^) in Fig. 8. Tunneling processes are expected to 
occur by the coherent resonant tunneling process in the solid line region. On the broken lines, the 
contribution from the processes with the inelastic spin excitation will be not small. The conductance 
show relatively large difference from Gf(<^) in this region. 
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Fig. 8. The Kondo temperature rK(v3) as a function of the magnetic flux 79 at ed = —0.025. Inset data are the 
magnetic excitation spectra at zero temperature. The energies oj of the peak positions of the spectra are used for 
the definition of TK{<fi)- Tk(v ~ 0) is extremely low because the hybridization intensity of the odd-channel, Ao, is 
proportional to sin^((^/4) and is very small for i/j ~ 0. Tk{(p) increases as (p increases from (p — to tv. 



When the temperature T rises beyond T7, the conductance near (^/vr ~ 1 gradually increases and 
the AB oscillations become relatively smaller. When the temperature increases further through 
Tio = 1.1 X 10^^, the conductance near ip/ir ^ gradually decreases and the AB oscillations 
become very small. The spin states of both even and odd channels gradually change from non- 
magnetic to magnetic state for (/j/vr ~ 1 as temperature rises up. 

We can conclude that the behavior of the conductance strongly reflects the spin state in the 
quantum dots. At this place we comments on the difference of G{(p) at very low temperature 
T = Ti from the curve of Gf(v') at exactly zero temperature. The conductance G^if = 0) is less 
than 0.1 X 2e^//i, but G{ip = 0) is almost 2e'^/h even at very low temperature T = Ti. In our 
calculation, it is very difficult to get the results which agree with Gf(v3) in the region yj/vr < 0.15, 
because we must continue the numerical computation to extremely low energy region. We have 
checked that calculations up to low energy region lo ^ gives consistent result 

with Gf(</?) for ip/iT = 0.117. (See Fig. 9). But it is impossible to continue such calculation to 
ip/iT < 0.117 cases. In the ip/ii < 0.15 cases at T = Ti, the contribution from the coherent resonant 
tunneling process of the odd-channel will be small because Tk(93) is less than Ti. In such cases, 
the conductance will have the value near 2e^ jh as seen from Fig. 3. We think that abrupt increase 
of G(}fi ~ 0) at T = Ti from G^i^p ~ 0) really reflects the crossover of states due to the increase of 
the temperature. 
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Fig. 9. The excitation spectra p" {uj)/uj at the temperature T = Ti = 5.8 x 10"'' (left figure) and T = 2.1 x 10"* 
(right figure) at (p/n = 0.117, £d = —0.025. The spectrum P {io)/Lu at T — 2.1 x 10"** gives consistent result with 
Gf{(p/tv = 0.117) = 2.4 X 10"^ at £d = -0.025. 



Here, we note the meaning of the spectrum P {uj)/uj. This quantity is related to the fluctuation 
of the bias voltage. It shows peak structure at energy comparable to the width of the Kondo 
resonance at the Fermi energy. It shows two peaks structure when (/? is small, as clearly seen in 
Fig. 9. These peaks seem to relate to the Kondo resonance in each channel. 
3.2.2 Shallow case (ea = O.OOO; 

Numerical results of the conductance G((/?) at = 0.000 for various T are shown in Fig. 10. 
The conductance G{(p) does not show rapid changes in this case, contrasted to the results shown 
in Fig. 5. We show the Kondo temperature Tk((/') in the = 0.000 case in Fig. 11. It is much 
higher than that of = —0.025 case in Fig. 8 at ~ 0. In the (?t.o,o) ~ 0.0 case, spin freedom 
on the odd-orbit will disappear and so the spin excitation mainly comes from the electrons in the 
even-orbit. Then the Kondo temperature Tk(v) is larger in this case. Even when the occupation 
number on odd-orbit increases due to the increase of the characteristic magnetic excitation 
energy is not so small because the electrons in both channels are in the mixed valence regime. So 
the variation of Tk((/') for change of is much smaller than that in the = —0.025 case. Therefore 
the conductance in Fig. 10 does not show the rapid change caused by the crossover of the spin state. 
In the temperature range T < Tg, the conductance shows simple cos 93 like behavior contrasted to 
the complicated if dependence of = —0.025 case. When the temperature T rises up than Tiq, 
the AB oscillations become smaller because the temperature becomes comparable to T]^{ip). 
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Fig. 10. The conductance G{ip) for various temperatures at Ed = 0.000. AB oscillations become smaller as temper- 
ature rises up than T ^ Tg. The phenomenon that G{(fi) rapidly changes with flux at low temperatures can not be 
seen in this case. 



§4. Summary and Discussion 

In this paper, we have investigated the temperature and the magnetic flux dependences of the 
conductance G for the AB circuit including two quantum dots. This quantity is expected to reflect 
the coherency of the electronic state, and is closely related to the 'Kondo effect in quantum dots'. 
The model Hamiltonian had the form of the two-channel Anderson model with even and odd 
orbits. The conductance was calculated by using the NRG method based on the formula without 
assumptions on tunneling processes. 

There were interaction terms between electrons in even and odd orbits, but the Kondo effects 
in the even-channel and the odd-channel occurred as if they were independent. When the dots 
system were in the Kondo regime near (/9 ~ 0, i.e., the occupation number in the odd-orbit satisfies 
('^o.o) ~ 1-0, the crossover of the spin state in the odd-channel occurred at a critical value of <p at low 
temperature. This is because the hybridization strength of odd-channel increases from the small 
value near 99 ~ to large value with increasing ip. AB oscillations had strong higher harmonics 
components reflecting the crossover of the spin state of the system. This result was contrasted 
to the single cos if dependence derived by Akera.El-' When the temperature rose up very high, the 
amplitude of the AB oscillations became smaller as expected. When the dots system were in the 
mixed valence regime near 99 ~ 0, i.e., (no,o) ~ 0.0, the drastic change of the conductance as a 
function of </? did not occur. 

At this place, we remark limitation of the assumptions used to derive the simplified expressions, 
eqs. (2.5), (2.6) and (2.7). The even and odd channels were separated out, and the hybridization 
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Fig. 11. The Kondo temperature TKif) as a function of the magnetic flux <^ at £d = 0.000 case. Inset data are 
the magnetic excitation spectra at zero temperature. The energies uj of the peak positions of the spectra are used 
for the definition of TKi^p)- The variation of T-K(tp) for the change of ip is smaUer when compared with result for 
Ed = —0.025 case. 



of odd-channel is proportional to sin((/7/4). This led the drastic change of the electronic state of 
the odd-channel near (/9 ~ 0. However, if the dot 1 and dot 2 are asymmetric under the interchange 
of them, electrons in both channels mix with each other. Therefore, the (/9-dependence near (/9 ~ 
will be blurred in the asymmetric case. The experimental work of Yacoby et al. has been done 
in this situation.llll© Extension of our calculation to the asymmetric cases will be given in near 
future. We stress that the present work treated the case that most drastic changes of the tunneling 
conductance are expected. 

Anyway, the conductance will show sharp (/9-dependence reflecting the crossover of the spin state 
of dots, when the Kondo effect is caused by the suitable variation of the dot potential energy. 
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Appendix: Formulation of the conductance in the Unear response theory 

We define an electric current as follows; 

J.-e-"^^' + <^°'. (A.l) 

where — e is the charge of an electron. The quantity A^a = ^[^;^a] is the time derivative of the 
electron number operator in the lead A, Na = Sfeo- ^lo-^fco-, and Nb is that in the lead B. 

We consider the current when the external perturbation term H = —NaVa — AbVb is added to 
H which is given by eq. (2.1). The expectation values, {Na) and (A'^b) are given within the linear 
response theory,!!^ 

{Na) = (taaVa + (JabVb, (A-2) 
{Nb) = (TbaVa + (TbbVb, (A-3) 

with 

afj^u = lim afj^uioj) (A-4) 
>o 

^ lim U r die-'' 

h [Jo 

xTr (^e-'f*'[iV,,/Jeq]e'f *'a^)| , (A-5) 

where 5 = 0+. Using properties of the trace operator and integration by parts, we get the following 
expression; 

a^,{io) = - {K^,{uj) - K^,{0)) , (A-6) 

lUJ 



with 



K^,{u:)^--l e'''+'-\[N,,N^{t)])dt. (A-7) 







The quantity K^^{lli) can be expressed by the following form; 

where K'^^{uj) and K'^^{uj) are given as follows by using the eigenstates of H (i.e., H\m) = Em\m)), 



z 



Em — En + flUJ 



[Em -En + hojf + {h5Y 

x{n\Ny\m){m\N^\n) (A-9) 
<.(-^^), (A-10) 
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n,m 

-ns 



X {n\Ni^\m){m\N^\n) 



The quantity u^i, is given by 



lim 



(A-ll) 
(A-12) 



(A. 13) 



We give the value — Va 
as foUows; 



-eV and — Vb = eV for external field, and obtain the electric current 



I = -r {o'AA + CBB — CTAB — Cba) " 2^. 
4 



The conductance G is given by 



2V 



2e2 P"(lo) 

— hm — . 

h u}-^o nuj 



(A-14) 

(A-15) 
(A- 16) 



with 



4 Z 



(n|iVA-iVB|m) 



X(5 {TUJ- {En - Em)) . 



(A-17) 



We note that the operator -/V^(jLt = A, B) can be expressed by the localized operators near the 
dots. For example, A'a — A^B is given as follows; 



2i 



iVA - iVe = - 5] y ( cos j4,a«a,fea 



(A- 18) 



In this place the operator J^ka'^iMi^ ~ S'^)> which is given by linear combination of all /c-state, 
is proportional to the localized orbit near the dots. Therefore the expression of conductance by 
using the quantity P" {oS) is suitable for the calculation based on the NRG method, in which the 
approximate eigenstates are obtained successively starting from the localized orbits. 
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